Asymptotic Regimes of Magnetic Bianchi Cos- 
mologies 
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Abstract. We consider the asymptotic dynamics of the Einstein- 
Maxwell field equations for the class of non-tilted Bianchi cosmolo- 
gies with a barotropic perfect fluid and a pure homogeneous source- 
free magnetic field, with emphasis on models of Bianchi type VIIq, 
which have not been previously studied. Using the orthonormal 
frame formalism and Hubble-normalized variables, we show that, 
as is the case for the previously studied class A magnetic Bianchi 
models, the magnetic Bianchi VIIq cosmologies also exhibit an os- 
cillatory approach to the initial singularity. However, in contrast 
to the other magnetic Bianchi models, we rigorously establish that 
typical magnetic Bianchi VIIo cosmologies exhibit the phenomena of 
asymptotic self-similarity breaking and Weyl curvature dominance 
in the late-time regime. 

Key words. Non-tilted magnetic Bianchi cosmologies. 

1 Introduction 

The influence of an intergalactic magnetic field on cosmological models has 
been investigated for over four decades both from a theoretical and observa- 
tional point of view. Cosmologists speculate that such a field could be pri- 
mordial in origin, that is, one that came into existence at the Planck time. 
Observational techniques rely on studying processes such as the tempera- 
ture distribution of the cosmic microwave background radiation (CMBR), 
primeval nucleosynthesis and the Faraday rotation of linearly polarized ra- 
diation emitted from extragalactic radio sources. Barrow et al [J derive an 
upper bound of < 3.4 x 10^^ (Oq/iso)^^^ gauss on the present strength of 
any spatially homogeneous primordial magnetic field^ based on data from 
the COBE satellite (fio is the present value of the density parameter and 
/i50 is the Hubble constant in units of 50kms~^ Mpc""'^). All observations 
to date only place an upper bound on the strength of such a magnetic field 
and hence are inconclusive as regards its existence. 

Any cosmological model which contains a magnetic field is necessar- 
ily anisotropic, since isotropy is violated by the preferred direction of the 
magnetic field vector. Consequently, one must analyze the Einstein field 
equations in models more general than the homogeneous and isotropic 
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■^For comparison, the strength of the Earth's magnetic field at the surface is approx- 
imately 0.5 gauss. 
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Friedmann-Lemaitre (FL) models. The simplest family of cosmological 
models that can admit a magnetic field are the so-called Bianchi cosmolo- 
gies, that is, models which admit a three-parameter group of isometries 
acting orthogonally transitively on spacelike hypersurfaces. The models 
are thus spatially homogeneous, but, in general, anisotropic. We assume 
that the models contain a barotropic perfect fluid whose four-velocity is 
orthogonal to the group orbits, and that observers comoving with the fluid 
measure a pure source-free magnetic field. We also assume that the per- 
fect fluid satisfies an equation of state p = (7 — 1)//, where 7 is constant 
and satisfies | < 7 < 2, the cases 7 = 1 (dust) and 7 = | (radiation) 
being of primary interest. We shall refer to solutions of the combined 
Einstein-Maxwell field euqations that satisfy the above properties as mag- 
netic Bianchi cosmologies. It is known that the field equations lead to 
restrictions on the Bianchi-Behr type of the isometry group, namely, that 
it is of types I, II, VIq or VIIq (in class A) or type III (in class B)^. 

Significant progress has been made in the study of magnetic Bianchi 
cosmologies. Collins [2| was the first to use techniques from dynamical 
systems theory to obtain qualitative results concerning the evolution of 
axisymmetric Bianchi I models under the assumption that the magnetic 
field is aligned along a shear eigenvector. More recently, LeBlanc et al |H| 
gave a qualitative analysis of the dynamics of magnetic Bianchi cosmologies 
of type VIq in their asymptotic regimes, that is, near the initial singularity 
and at late times. This work made use of the orthonormal frame formalism 
of Ellis and MacCallum ^ and Hubble- normalized variables (see ^H] , ch. 5 
and 6) . This formalism was also used in |H] and [Jj to give a similar analysis 
of magnetic Bianchi universes of types I and II. Most recently, Crowe |H] 
extended the results to magnetic Bianchi models of type III. There remains 
one class which has not been previously analyzed, namely magnetic Bianchi 
cosmologies of type VIIq. 

Our goal in this paper is to fill this gap by giving a qualitative analysis of 
the dynamics of magnetic Bianchi cosmologies of type VIIq in their asymp- 
totic regimes. Bianchi cosmologies of type VIIq are of interest because they 
represent anisotropic generalizations of the flat FL models. The asymptotic 
dynamics of the non-magnetic models at late times has only been analyzed 
in detail relatively recently (see j2()l ). It is worth comparing non-magnetic 
Bianchi cosmologies of group type VIIq with their counterparts of group 
type I, which are also anisotropic generalizations of the flat FL models. The 
non-magnetic Bianchi type I cosmologies are asymptotically self-similar at 
late times, that is, they are approximated by a self-similar solution at late 
times. This self-similar solution is in fact the flat FL solution, which means 
that the Bianchi I cosmologies undergo asymptotic isotropization. In con- 

•^We refer to Ellis and MacCallum |4| for this terminology. 
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trast, for values of the equation of state parameter 7 satisfying | < 7 < 2, 
the Bianchi VIIq cosmologies are not asymptotically self-similar at late 
times. Nevertheless, for values of 7 satisying 1 ^ 7 ^ |, they undergo a 
subtle form of isotropization: the rate of expansion isotropizes, but the in- 
trinsic gravitational field, as described by the Weyl curvature tensor, does 
not. This phenomenon has been referred to as Weyl curvature dominance 
(see [20]). One of our specific goals in this paper is to determine what effect 
a cosmic magnetic field has on the above-mentioned isotropization. The 
method that we use is a generalization of the analysis of the non-magnetic 
Bianchi VIIq and VIII models given in [201 ^^'^ 01 respectively. 

The plan of paper is as follows. In section |21 we present the evolu- 
tion equations for the magnetic Bianchi cosmologies of type VIIq using the 
orthonormal- frame formalism and Hubble- normalized variables. Section |21 
contains the main result concerning the dynamics in the late-time regime, 
namely theorem l^.ll and corollarv l3.ll which give the limits of the Hubble- 
normalized variables and of certain physical dimensionless scalars, thereby 
describing the asymptotic dynamics at late times. In section 01 we exam- 
ine the singular asymptotic regime and show that typical models exhibit 
an oscillatory singularity. We conclude in section with a discussion of 
the cosmological implications of our results and give an overview of the 
asymptotic dynamics of the magnetic Bianchi cosmologies, noting that the 
present paper completes the picture. 

There are three appendices. Appendix A contains the proof of the fact 
that magnetic Bianchi VIIq universes are not asymptotically self-similar 
at late times. Appendix B fills in some of the technical details of the 
proof of theorem 13.11 Finally, in appendix C we give expressions for a 
dimensionless scalar formed from the Weyl curvature tensor in terms of 
the Hubble-normalized variables. 

2 Evolution equations 

In this section we give the evolution equations for magnetic Bianchi cos- 
mologies of type VIIo. As described in |H] (pg. 517), we use Hubble- 
normalized variables 

(E+,E_,7V2,iV3,?^), (2.1) 

defined relative to a group-invariant orthonormal frame {fia}, with Bq = u, 
the fluid 4-velocity, which is normal to the group orbits. 

The variables E± describe the shear of the fluid congruence, the A^2,3 are 
spatial connection variables which describe the intrinsic curvature of the 
group orbits and the variable Ti describes the magnetic degree of freedom. 
The magnetic Bianchi VIIq cosmologies are described by the inequalities 
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H > and N2N3 > 0. Without loss of generality, we assume 

^2 > 0, 7V3 > 0. (2.2) 

It is convenient to define 

N+ = i(7V2 + ^3), = ^{N2- iVa), (2.3) 

and replace (|2.1|) by the state vector 

(I]+,I]_,7V+,iV_,7i). (2.4) 

The restrictions (12.21) become 

7V+ > 0, Nl- ml > 0, H > 0. (2.5) 

The state variables H2.1II and (|2.4|l are dimensionless, having been nor- 
malized with the Hubble scalar'' H, which is related to the overall length 
scale £ by 

H=j, (2.6) 

where the overdot denotes differentiation with respect to clock time along 
the fundamental congruence. The state variables depend on a dimensionless 
time variable r that is related to the length scale i by 

e = ioe\ (2.7) 

where io is a constant. The dimensionless time r is related to the clock 
time t by 

as follows from equations H2.6|l and (|2.7|l . In formulating the evolution 
equations we require the deceleration parameter q, defined by 

(2.9) 

and the density parameter fJ, defined by 

= -^. (2.10) 

*On account of 12. 6i . H is related to the rate of volume expansion 9 of the funda- 
mental congruence according to H = ^6. We note that all variables of LeBlanc et al |S] 
are normalized with 8. 
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fih = ^{hl + hl + hl), (2.12) 



We also find it convenient to introduce the magnetic density parameter £7^, 
defined analogously by 

ft. ^ ^. (2^11) 

where Hh is the energy density of the magnetic field. We note that fj.h is 
given by 

where the ha, a = 1, 2, 3, are the components of the magnetic field intensity 
relative to the spatial orthonormal frame {gq}, which has been chosen so 
that 

ha ^ {hi, 0,0) (2.13) 

(see El). 

A complete derivation of the evolution equations for the variables H2.4|l . 
which arise from the combined Einstein-Maxwell field equations, is provided 
in [HI (see section 2). These evolution equations read^ 

= (g - 2)2+ - 2N^ + 
E'_ = (g-2)I]_ -2N+N-, 

A^l = (g + 2I]+)iV+ + 6S_iV_, (2.14) 
N'_ = iq + 2S+)iV_ + 2S_iV+, 
n' = (g-2E+- 1)H, 

where 

q = 2(j:l + j:'i) + lH^ + ^{3j-2)n, (2.15) 

fl = l-T.l-j:^_- Nl- iW^, (2.16) 

and ' denotes differentiation with respect to r. For future reference we also 
note the evolution equation for fl: 

n' = [2q - - 2)]n, (2.17) 

and the expression for the magnetic density parameter 

= \n\ (2.18) 

in terms of the Hubble- normalized magnetic field intensity 7i — hi/H, 
which follows from l(TTT|l . and itTT^ . 

^ These evolution equations are essentially the same as those given in |H| for magnetic 
Bianchi VIq models, apart from a numerical factor multiplying Ti^. The difference 
between Bianchi VIIq and Bianchi VIq models lies in the restrictions that define the 
state space: the quantity — 3N^ is negative for Bianchi VIq models, in contrast to 
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The physical requirement f2 ^ in conjunction with H2.3I) imphes that 
the variables S± , N- and H are bounded, but places no restriction on 
itself. In fact, it will be shown in appendix A (see proposition IA.1|I that if 
> and I < 7 < 2, then for any initial conditions 

lim N+ = +00. (2.19) 

r — *-j-oc 

The first step in analyzing the dynamics at late times (r — > +oo) is to in- 
troduce new variables which are bounded at late times and which enable us 
to isolate the oscillatory behaviour associated with S_ and N-. Motivated 
by [20], we define 

I]_=i?cosV', N-^RsiniP, M = (2.20) 

where i? ^ 0. 

In terms of the new variables (E+, i?, 7i, M, V'), the evolution equa- 
tions (|2.14l) have the following form 

= (0 - 2)E+ -R^ + + (1 + S+)i?2 cos 2^), 
i?' = [Q + - I + {R^ - I - cos 2if\ R, 
7i:'= [Q-2S+-l + i?2cos2i/']H, (2 21) 

M' = - [g + 2S+ + i?2(cos2i/' + 3Msin2V')] M, 

= ^[2+(l + E+)Msin2V'], 



where 
and 



^ 2Y.\ + R^ + }re + - 2)VL, (2.22) 



n = l-i:\- R^ -VLh. (2.23) 
The evolution equation for 17 becomes 

n' = [2Q - (37 - 2) + 2i?2 cos 2?/'] ^. (2.24) 

The restrictions (12. 5() are equivalent to 

3M2i?2 sin2 ^ < 1, Af > 0, i? ^ 0, H > 0. (2.25) 

3 Limits at late times 

In this section we present a theorem which gives the limiting behaviour as 
r +00 of the magnetic Bianchi VIIo cosmologies when the equation of 
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state parameter 7 satisfies | < 7 < 2. As a corollary of the theorem, we 
obtain the limiting behaviour of certain dimensionless scalars that describe 
physical properties of the models, namely the density parameter Q, defined 
by (|2.1()|l . the magnetic density parameter H.^, defined by (|2.11(l . the shear 
parameter S, defined by 

- (3-1) 

where cjah is the rate-of-shear tensor of the fluid congruence, and the Weyl 
curvature parameter W, defined by 

= ^"^"g+f"''" , (3.2) 

where E^h and H^h are the electric and magnetic parts of the Weyl tensor, 
respectively (see JHl, Pg- 19), relative to the fluid congruence. 

In terms of the Hubble-normalized variables, the shear parameter is 
given by 

I]2 = + i?2 cos^ V-, (3.3) 

which follows from H2.20|l in conjunction with equation (6.13) in The 
formula for the Weyl curvature parameter is more complicated and is pro- 
vided in appendix C. 

The main result concerning the limits of S+, i?, H, and M is contained 
in the following theorem. Some of the results depend on requiring that the 
model is not locally rotationally symmetric^ (LRS). 

Theorem 3.1. For all magnetic Bianchi cosmologies of type VIIq that 
are not LRS and with density parameter 17 satisfying > 0, the Hubble- 
normalized state variables (I]+,i?, 7i,M) satisfy^ 

if |<7<|, 
S+,i?,H,Af) = <! f0,A/ifc,V2A;,0) , if 




liin _(S+,i?,H,Af) = <j (0, ^|fc,%/2A:,0j , if 7 = |, (3.4) 

^,%/2,0), if |<7<2, 



and 




(3.5) 



^See, for example, 1191 . pg. 22. We note that the LRS magnetic Bianchi VIIo models 
are described by the invariant subset S_ = N— = 0, equivalently, R = 0. Since LRS 
models of Bianchi type VIIq also admit a group G3 of isometrics of Bianchi type I, we 
do not consider them in detail here. 

^The limits in the case 7=5 were conjectured by Sam Lisi. We thank him for 
helpful discussions. 
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where fc G (0, 1) and L > are constants that depend on the initial condi- 
tions. 

Proof. It foUows immediately from H2.19|l and H2.20|l that 

lim M = 0. (3.6) 

T — *-t-00 

Furthermore, since S+ is bounded, it follows fr'om the ip evolution equation 
in that 

lim "0 = +00. (3-7) 

r — >+oo 

The trigonometric functions in the DE H2.21|l thus oscillate increasing 
rapidly as r — > +00. In order to control these oscillations, we introduce 
new gravitational variables S+, R and according to^ 

S+ = S+ - i(l + i:+)R^M sin2V', 
R = R[l-\{R^ -l-Y.+)Msm 2^)] , (3.8) 
7Y = 7i:[l-ii?2Msin2?/'] . 

These new variables are defined so as to 'suppress' the rapidly oscillating 
terms which may not tend to zero as r — > +00. The evolution equations 
for these "barred" variables, which can be derived from (|2.21() and H3.8|l . 
have the following form 

= -(2 - Q)S+ - i?2 + 17^2 + AfBs+, 
i?' = (g + S+ - 1 + MBj{)R, (3.9) 
7Y' = (Q - 2E+ - 1 + MBfi)n, 

where 

Q^2tl + R^ + + i(37 - 2) (1 - - i?2 „ 17^2) ^ (3.10) 

and the B terms are bounded functions in E+, R, Ti and in M and tp for 
r sufficiently large. The essential idea is to regard M and ■0 as arbitrary 
functions of r subject only to H3.6|) . Thus, (|3.9|) is a non-autonomous DE 
for 

S= (S+,i?,7i), 

of the form 

x' = fix)+g{x,T), (3.11) 

where 

gix,T) = MiT)iB^^,RBj,,nB^), (3.12) 

*We are motivated by the analysis in (S| (see equations (3.8)) and I2UI (see equa- 
tions (B.4)). 
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Table 1: Limits of the Hubble- normalized scalars fl, fit, S and W at late 
times for magnetic Bianchi cosmologies of type VIIq. 



Range of 7 


n 




E2 




w 


i<7<l 


1 













7=1 


1 










L / 


1<7<| 


1 










+00 


7=1 






(0, 




+00 


|<7<2 





1 
3 


(0, 


i) 


+00 



tThe components in the parentheses are the liminf and Umsup. The parameter k is 
the parameter that appears in theorem l3.1l 

and f{x) can be read off from the right-hand side of (|3.9|) . Since 

lim g{x, t) = 0, 

r — ^-t-00 

as follows from H3.6|l . the DE H3.9|l is asymptotically autonomous (see |12)). 
The corresponding autonomous DE is 

& = (3.13) 

where 

X = it+,R,'H). 

Using standard methods from the theory of dynamical systems, we first 
show that the limits of the "hatted" variables correspond to those limits 
stated in the theorem. Details are provided in appendix B.l. We then 
use a theorem from J^l (see theorem IB. II in appendix B) to infer that the 
solutions of the non- autonomous DE (|3.11|) have the same limits as the so- 
lutions of the autonomous DE H3.13|l . Details are provided in appendix B.2. 
The hmit of X — (S]_|_, i?, 7Y) follows immediately from this result in con- 
junction with the definitions (|3.8() . Finally, the limit l|3.5ll concerning the 
ratio M/R is obvious when | ^ 7 < 2, since R ^ 0. The more complicated 
case when | < 7 < | is treated in appendix B.3. □ 

Corollary 3.1. The limits as r ^ +00 of the density parameter il, the 
magnetic density parameter flh, the shear scalar E and the Weyl curvature 
scalar W, for all magnetic Bianchi cosmologies of type VIIo that are not 
LRS and with satisfying il > 0, are as given in tabled 

Proof. These results follow directly from theorem |3.1l and equations H2.23II . 
(ESI, lfCl3)l and (|Cl4)l . Moreover, if 1 < 7 < 2, it follows from l|Cl3)l and 
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ljC.4|) that since T,^, R and Ti. arc bounded, and lim^^+oo M/R = 0, that 



W^—[l + 0{M)], 

as r ^ +00. □ 
We conclude this section by discussing the physical interpretations of 
theorem 13.11 and its corollary. Like the non-magnetic Bianchi VIIq mod- 
els, the magnetic Bianchi VIIo models are not asymptotically self-similar 
at late times, since the orbits in the Hubble-normalized state space do 
not approach an equilibrium point of the evolution equations. This phe- 
nomenon is accompanied by Weyl curvature dominance, characterized by 
the divergence of the Weyl curvature scalar W which describes the intrinsic 
anisotropy of the gravitational field. For non-LRS models with 1 < 7 < 2, 
W is unbounded as r ^ -l-oo. 

The shear scalar S quantifies the anisotropy in the expansion of a cos- 
mological model. We see that for | < 7 < |, the models isotropize at 
late times in the sense that limT-^+oo S = 0, as is the case for the corre- 
sponding non-magnetic Bianchi VIIq models (see ^01, theorem 2.3). The 
key difference between the magnetic and non-magnetic models occurs when 
the matter content is a radiation fluid (7 — |): the presence of a magnetic 
field prevents shear isotropization, in the sense that £ does not tend to 
zero at late times. 



4 The singular asymptotic regime 

In this section we show, by combining numerical experiments with ana- 
lytical considerations, that generic non-LRS magnetic Bianchi VIIq cos- 
mologies with I < 7 < 2 exhibit an oscillatory approach to the initial 
singularity, as do the other class A magnetic Bianchi models. 

As is the case for the previously studied class A magnetic Bianchi mod- 
els, the behaviour into the past for the magnetic Bianchi VIIq models is 
necessarily complicated since none of the equilibrium points of the evolution 
equations H2.14|l are local sources. It is well-known that in the dynamical 
systems approach, the Kasner circle /C plays a primary role in determining 
the dynamics towards the singularity, since its local stability enables one 
to predict whether the singularity in a given class of models is oscillatory^. 
For the present class of models, /C is the set of equilibrium points described 
bylO 

S+ + = 1, = Ns^n^o. 

^In a recent paper 1181 . it has been shown that the Kasner circle also plays this role 
in cosmological models without symmetry. 

^"in discussing the singular asymptotic regime, it is more convenient to use the spatial 
connection variables A^2 and N3, rather than and A''— (see equation 12.31 '). 
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Figure 1: The arrays — h +, etc. give the signs of the eigenvalues A-h, 
and Aatj in that order. The variables listed next to each of the three arcs 
indicates which of the variables TC, is growing into the past. 

A local stability analysis shows that the Kasner equilibrium points are 
saddles in the Hubble-normalized state space. Apart from three exceptional 
points (labeled Ti, T2 and in figure^, the equilibrium points of IC have 
a one-dimensional unstable manifold into the past. Figure ^ shows the 
signs of the eigenvalues on JC associated with the variables H, N2 and 
N3 (a negative eigenvalue indicates instability into the past) and which of 
these three variables are increasing into the past in a neighbourhood of the 
Kasner circle. 

It turns out each unstable manifold on IC is asymptotic to another 
Kasner point. In other words, the unstable manifold is a heteroclinic orbit 
of /C, i.e. an orbit which joins two Kasner points. These unstable manifolds 
provide a mechanism for a cosmological model to make a transition from 
one (approximate) Kasner state to another, as it evolves into the past. 
Figure |2] shows the projections in the S+S_-plane, of these families of 
heteroclinic orbits, that join two Kasner points. The families are described 
as follows: 

Sn- 1^1 + 1:1 + Ire = 1, 7^ > 0, iVs = A^a = 0, 

Sn, : + + j^A^I = 1, 7V2 > 0, Afg = W = 0, (4.1) 

Sn, ■■ ^1 + ^1 + ^Nl = 1, A^a > 0, iVs - H = 0. 

The heteroclinic orbits on 3 describe the familiar vacuum Bianchi II 
Taub models, while the orbits on Sn describe the Rosen magneto-vacuum 
models (see jH|, pg. 531). 

Numerical experiments suggest that for generic non-LRS models, after 
an initial transient stage, the orbit approaches a point on IC. The direc- 
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Figure 2: The projections of the Rosen and Taub orbits joining points on 
the Kasner circle /C. The arrows show evolution into the past. 



tion of departure of the orbit is determined by the unique Rosen or Taub 
orbit through that point whereupon it shadows (i.e. is approximated by) 
this orbit until is approaches another point on /C and the process repeats 
indefinitely. In physical terms, the corresponding cosmological model is 
approximated by an infinite sequence of Kasner vacuum models as the sin- 
gularity is approached into the past, the so-called Mixmaster oscillatory 
singularity. This behaviour motivates the following conjecture concerning 
the past attractor . 

Conjecture 4.1. The past attractor is the two-dimensional invariant set 
consisting of all orbits in the invariant sets S-^, Sjq^ and Sjq.^ (see figure E)) 
and the Kasner equilibrium points, i.e. 

A^ =Sh^ Sn2 U Sn^ U JC. (4.2) 



This conjecture can be formulated in terms of limits of the state vari- 
ables as follows. Referring to (j2.16(l and (j4.1|l . we see that the set A" is 
defined by = and 

iVaH = NsH = N2N3 = 0. 

It follows from monotone function arguments (see the comment at the end 
of appendix A) that 

lim N2N3 = 0, 

r — ^ — 00 

and, moreover, that and N3 are bounded in the singular regime. Thus 
our conjecture concerning the past attractor can be formulated as 

lim n^O, lim = 0, lim N^H = 0. 

r — >~oo T — * — 00 r — > — QO 

Note that for a generic orbit, limr^-ooiTi-, N2, N3) does not exist. 
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Table 2: Shear, spatial curvature and magnetic degrees of freedom in 
class A Bianchi cosmologies. 



Bianchi type Shear Spatial curvature Magnetic field 

I 2 3 

II 2 1 2 
VIo, VIIo 2 2 1 
VIII, IX 2 3 



5 Discussion 

With the appearance of the present paper there is now available a com- 
plete description of the dynamics of magnetic Bianchi cosmologies^^ with a 
perfect-fluid matter content, in the two asymptotic regimes. We now give 
an overview of the properties of these models, in order to highlight the 
role of a primordial magnetic field in spatially homogeneous cosmological 
dynamics. The possible Bianchi types and relevant references are given in 
the introduction. We emphasize the so-called class A models (in the ter- 
minology of Ellis and MacCallum that is, those of Bianchi types I, II, 
VIo and VIIq. For each of these types the Hubble- normalized state space 
is five-dimensional, but, as indicated in table |21 they differ as regards the 
number of degrees of freedom associated with spatial curvature and with 
the magnetic field. In this table we have also listed Bianchi types VIII and 
IX, which do not admit a magnetic field, for comparison purposes. 

All models in table El display an oscillatory approach to the singularity, 
described by a two-dimensional attractor in the Hubble-normalized state 
space, familiar from the non-magnetic Bianchi VIII and IX models (see |19j . 
pp. 143-7). The essential point is that the magnetic field mimics spatial 
curvature in that it destabilizes the Kasner circle of equilibrium points. 
Note that the sum of the number of spatial curvature and magnetic degrees 
of freedom is three, in all cases in table [3 

As regards the late-time dynamics of the magnetic cosmologies, there 
is a fundamental difference between the Bianchi VIIq models considered in 
the present paper and the Bianchi I, II and VIq models considered earlier 
(0) |Z|)j as follows. The Bianchi I, II and VIo models are asymptoti- 
cally self-similar, in the sense that each model is approximated by an exact 
self-similar solution, while the models of type VIIq are not asymptotically 

^^We emphasize that we are restricting our attention to Bianchi cosmologies that are 
non-tilted, in the sense that the fluid four-velocity is orthogonal to the group orbits. 
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Table 3: Limiting values of Q, Q^, S and W as t +00 for magnetic 
Bianchi cosmologies with a dust fluid (7 = 1). 



Bianchi type Q 



I 1 

II ± ^ 

16 64 2048 

VIot 1(1 -fc') Ik' ^ ^(l + 2fc2)(2 + fc^) 

VIIo 1 i > 



tThe parameter k satisfies < fc < 1 and depends on the initial conditions. 



self-similar. This difference is essentially a consequence of the fact that the 
Hubble-normalized state space of the Bianchi VIIq models is unbounded. 
Another feature of the magnetic cosmologies is that the asymptotic dynam- 
ics at late times depends significantly on the equation of state parameter 
7. We illustrate this dependence in tablesOand^where we give the limits 
at late times of il, fth, S and W for the two physically important cases, 
dust (7=1) and radiation (7 = |). It is worthy of note that if 7 < |, the 
Bianchi I models isotropize in all respects (S ^ 0, W ^ and — > 0) 
while the Bianchi VIIq models isotropize as regards the shear and the mag- 
netic field (E ^ and nh^O). 

A cosmic magnetic field also affects the local stability of the equilibrium 
point that corresponds to the flat FL solution^^. For non-magnetic models, 
the flat FL equilibrium point is typically a saddle point, having both a 
non-trivial stable manifold and a non-trivial unstable manifold. The shear 
degrees of freedom generate the stable manifold while the spatial curvature 
degrees of freedom generate the unstable manifold. The stable manifold 
leads to the phenomenon of intermediate isotropization, i.e. a model can 
evolve to become arbitrarily close to isotropy over a finite interval of time. 
The unstable manifold leads to models with an isotropic singularity, i.e. 
models which are highly isotropic near the initial singularity, but which 
subsequently develop anisotropies. The effect of a primordial magnetic 
field on these phenomena depends on the equation of state parameter 7. If 
7 < |, the magnetic field increases the dimension of the stable manifold, 
leaving the dimension of the unstable manifold unchanged, thus increasing 
the likelihood of intermediate isotropization. On the other hand, if 7 > |, 
the magnetic field increases the dimension of the unstable manifold by one, 
leading to magnetic models with an isotropic singularity. 

We conclude by giving some suggestions for future research. Firstly, it 

^^In the present paper, this oquiUbrium point is given by T,± = 0, N± = 0, W = 0. 
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Table 4: Limiting values of fl, flh, S and W as r ^ +00 for magnetic 
Bianchi cosmologies with a radiation fluid (7 = |). 



Bianchi type 


n 








I 


1 











II 


3 


1 


1 


21 


4 


12 


12 


144 


VIo 





3 

8 


1 

16 


81 
128 


VII0+ 








+00 



^^The parameter k satisfies < A; < 1 and depends on the initial conditions. The 
components in the parentheses for correspond to its hminf and hmsup. 

would be of interest to investigate the asymptotic dynamics of spatially in- 
homogeneous cosmological models in the presence of a primordial magnetic 
field, in order to determine which features of magnetic Bianchi cosmologies 
occur in models without symmetries. The recent paper |18| on Go cosmolo- 
gies would provide a suitable framework for such an investigation. One step 
has been taken in this direction by Weaver et al who investigated a 
family of inhomogeneous cosmologies that generalize the magnetic Bianchi 
VIq cosmologies, and provided numerical evidence that the singularity is 
oscillatory. 

Secondly, the work of Barrow et al 1^ referred to in the introduc- 
tion leads to an upper bound on the magnetic density parameter fi/i of 
order 10~^° for spatially homogeneous magnetic fields. It would be of 
interest to what extent this upper bound would be weakened within the 
class of spatially inhomogeneous magnetic cosmologies. The analysis of the 
anisotropics in the CMBR by Maartens et al |^ would probably provide a 
suitable framework for such an analysis. 

We finally comment briefly on other recent work on magnetic fields in 
cosmology, which has focused on the potential dynamical effects of a pri- 
mordial magnetic field in a perturbed FL cosmology (jJOIj d]; dli 
jl5| . jl7|'l. or in a perturbed Bianchi I cosmology (JHI)- This work, which 
makes use of the Ellis-Bruni covariant and gauge-invariant method for an- 
alyzing cosmological density perturbations, complements the results in our 
paper and related ones, which focus on the dynamics in the asymptotic 
regimes of magnetic Bianchi cosmologies, and on the likelihood that such 
a model will evolve to be close to FL. Extending the dynamical systems 
analysis of magnetic Bianchi cosmologies to spatially inhomogeneous mag- 
netic cosmologies may help to bridge the gap between these two bodies of 
work. 
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Appendix A 

In this appendix we prove H2.19|l . concerning the limit of the Hubble- 
normalized variable iV+ in the latc-time regime. This result is restated 
as proposition lA. II below. 

Proposition A.l. For all magnetic Bianchi cosmologies of type VIIq that 
are not LRS"^'^, with equation of state parameter 7 subject to | < 7 < 2, 
and density parameter satisfying > 0, the Hubble-normalized variable 
A^-i- satisfies 

lim N+ = +00. (A.l) 

r — >+oc 



Proof. The proof is similar to the proof of the corresponding result for non- 
magnetic Bianchi VIIq cosmologies (see theorem 2.1 and equation (A. 5) in 
|20|). in that it makes use of monotone functions and the so-called mono- 
tonicity principle (see chapter 4 in ^Hl)- There are two cases depending on 
the value of 7. 

Case 1. I < 7 ^ 1 

As in [201, we consider the function 

with V — -1(37 — 2). The evolution equations (|2.14(l imply that 

Z[ _ + vf + [l - v'')Y?_] [1 + v){l - Av)n'^ 

Zi ^ 1 + VY.+ ^ 3(l + i;S+) 

Since < w ^ j in this case, it follows that Zi is monotone increasing 
and the result (|A.1|) follows as in the non-magnetic case (see appendix A 
in [2n|)- 

Case 2. 1 < 7 < 2 

^^The result also holds for the LRS models. The proof is similar to the non-LRS case; 
we omit the details in this paper. 



16 



We give a proof by contradiction. Suppose that (|A.1|I does not hold. Since 
the remaining variables are bounded, it follows that for any point x in the 
state space the w-limit set oj{x) is non-empty. 
Consider the function 



S2 



(A.3) 



which satisfies < Z2 < +00 on the invariant set S defined by 
N+ > 0, Nl- 3Nl > 0, T.^_+Nl>0, f7 ^ 0, H^O 

The evolution equations 12.1411 imply that 

_ -4(1 + S+)S^ 
Z2 " + iV2 ■ 

It follows that Z2 is decreasing^^ along orbits in S. We can now apply the 
monotonicity principle. By (A.3) the set S \ S (the set of boundary points 
of S that are not contained in S) is defined by one or both of the following 
equalities holding 

Nl - 3Nl =0, ^l+Nl ^ 0. (A.4) 

It now follows that for any x Cz S, the w-limit set uj{x) is contained in the 
subset of 5 \ 5 that satisfies lim^^g ^2(2/) 7^ +00, where s E S \ S and 
y G S. On account of (|A.2p and (|A.4|I we conclude that 



uj{x) c {a; I S- iV_ = 0}. (A.5) 
We can further restrict the possible w-limit sets by considering the function 

^ (iV2_ 3^2)^2' 

which satisfies 

= -6(7 -1)^3, 

as follows from 12.14(1 . We immediately conclude that limr^+00 Z3 = and 
hence that limr^+00 = 0. In conjunction with ((A. 5(1 . this result implies 
that a; (a;) C Si, where 

Si = {x\T.^ = = Q = 0}. 

The only potential w-limit sets in 5*1 are equilibrium points, since"*^^ lim^^+oc 
+00 for all other orbits in S'l. The equilibrium points are 

-'^■'No orbit in S satisfies S4. = — 1 for all r. 

^^On Si the evolution equation for A''+ reduces to N', = (1 + T:+)^N+. 
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i. E+ = 1, 7V+ = H = 0, (an isolated point) 

ii. = H = 0, iV+ > (aline). 

No orbit with H > 0, > and 1 < 7 < 2 can be future asymptotic to any 
of these equilibrium points, since fi' > in a neighbourhood of any of these 
points as follows from (|2.17|) and (|2.15|l . Thus we have a contradiction of 
the fact that uj{x) (j>, and as a result, ljA.l|) holds in case 2. □ 

Comment. The monotone function l|A.2p also provides useful information 
about the past asymptotics of magnetic Bianchi cosmologies. From the 
monotonicity principle, we can conclude that 

a{x) C {cc I iV^ - 3Nl = 0}, 

for any x d S. Therefore, in contrast to the late-time regime, A^_|_ is 
bounded towards the initial singularity. 

Appendix B 

In this appendix we fill in the details of the proof of theorem 13.11 The 
proof of this theorem relies on a result of ^21 (see corollary 3.3, pg. 180) 
concerning asymptotically autonomous DEs, stated as theorem IB . II below. 
Consider a non-autonomous DE 

x' ^ fix)+gix,T), (B.l) 

and the associated autonomous DE 

x' = fix), (B.2) 

where / : Z? -> M", g : D x M ^ K" and Z? is an open subset of R". It is 
assumed that 

Hi: lim g(w{T),T) — for every continuous function w : [ro,-|~oo) — > 

r — >+oo 

D 

and 

H2'. any solution of HB.1() with initial condition in D is bounded for r ^ 
To, for some tq sufficiently large. 

Theorem B.l. li Hi and H2 are satisfied and any solution of (|B.2p with 
initial condition in D satisfies 

lim x{t) ~ a, 

then any solution of IjB.ip with initial condition in D satisfies 

lim x{t) — a. 

We make use of this theorem in appendix B.2. 
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Appendix B.l 

We now deduce the limits at late times of , Rjli,). The components of 
the DE H3.13|l . x' = /(&), are given by 

R' ^{Q + ±+- 1)R, (B.3) 

n' = {Q- 2t+ - i)n, 



where 



Q = + i?2 + 17^2 _^ 1(37-2)17, (B.4) 
Cl = l-t.l-R^ - \re. (B.5) 



One can also form an auxilliary DE for using ifO)) and (P;5|) to find 
that 

f2' = [2Q- (37- 2)]f2. (B.6) 
We consider the state space S of the DE (|B.3p defined by the inequalities 



i?>0, n>Q, n>{). (B.7) 



These inequalities in conjunction with ljB.5|) imply that the state space S 
is the interior of one quarter of an ellipsoid. Understanding the dynamics 
on the two-dimensional invariant sets 5^, Sj^ and Sj^, the closure of their 
union defining the boundary of S*, will be crucial in our analysis. These 
sets are defined by the following restrictions: 





n 


= 0, 


i? > 0, 


7^ > 0, 


Sr-- 


R 


= 0, 


H>0, 


n>o, 


'■ 


H 


^ 0, 


i? > 0, 





The DE (|B.3|I admits a positive monotone function 

Z^^-^, (B.8) 
Rm^ 

which satisfies 

Z' = 3(4-37)Z (B.9) 

on the set S. Thus, if 7 7^ | there are no equilibrium points, periodic orbits 
and homoclinic orbits in S (see |19j . proposition 4.2). It is immediate upon 
integrating ljB.9p and using the boundedness of R and H that for any x G S 

ioix) C Sj^U Sfy, if |<7<|, (B.IO) 
uj{x) CS^, if I < 7 =s: 2. (B.ll) 
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Figure 3: Orbits in the invariant set 



We now consider the case | < 7 ^ 2. The flow on the invariant set S^^ 
is depicted in figureOl which shows the projection of the surface = onto 
the E+i?-plane. The essential features are the existence of three equihbrium 
points 

K±: (E+,7?,7i) = (±1,0,0), 
V: (E+,i?,7i)= (^0,y^,\/2^, 

in which he on the boundary of S^, and the fact that there are no 
periodic orbits on 5*^^. The latter can be established by the existence of a 
Dulac function A on given by 

A = R-^{1 - S+ - R^)~^ 

(see theorem 4.6, pg. 94). Thus, the only potential w-limit sets in 
S'q are the equilibrium points K^, V and the heteroclinic sequence (K~ — > 
K+ K~) and hence for any 5; € 5^, the w-limit set is one of these four 
candidates. The point K+ can be excluded since it is a local source in 
S. Moreover, the point K~ can be excluded by considering the evolution 
equation for 7i, which is of the form 

n' = h{t+,R,n)n. 

Since /i(K-) = /i(-l,0,0) = 3 and = at K", it follows that lim^^+oc 7^ 
and hence that an orbit in S cannot be future asymptotic to K^. This 
leaves the equilibrium point V and the heteroclinic sequence (K^ K+ 
K~) as the remaining candidates for the w-limit set in S^. The latter can 
be excluded since V is a local sink in S and hence lu{x) — V for any x E S^. 
On account of (|B.lip . we thus conclude that for any x £ S, 

lim (S+,i?,7i) = ( 0,./|,V2) , if |<7<2 (B.12) 
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n 




The case f < 7 < | can be treated in a similar fashion by analyzing 
the dynamics on the invariant sets Sj^ and S^. It follows that 

lim (E+,i?,7i;) = (0,0,0), if |<7^t (B-13) 

Finally, we consider the case 7 = | ■ We first observe that the VL evolu- 
tion equation (|B.6|) restricted to a radiation fluid reduces to 

It follows immediately from the LaSalle invariance principle (see the- 
orem 4.11, pg. 103) that 

uj{x) C {& I S+ = 0}, (B.14) 



for any x £ S. By (|B.9p the function Z defined in (jB.8|) describes a 
conserved quantity 

= k, (B.15) 
where > is a constant that depends on the initial condition. We see 



that for all A: > the surfaces described by l|B.15(l foliate the state space 
S and intersect the boundary O = at i? = and H = (see figure EJ. 
When 7=1, the DE (|B.3|I has a line L of equilibrium points given by 



(t+,R,n)= {Q,J^k,V2k] , fee (0,1) 



3 

It can show that for each fc > 0, the two-dimensional invariant set defined 



by (jB.15|l intersects the line L at precisely one point. Since this unique 
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point of intersection is the only equilibrium point on this invariant set 
which satisfies S+ = 0, it follows from the restriction l|B.14|l that any 
solution in S satisfies 

lim {t+,R,n)= (o,^/^k,V2k] , if 7 = |, (B.16) 
\ V J 

where fc G (0,1) is a constant which depends on the initial condition. 



Appendix B.2 

We now apply theorem lB . II using the results of appendix B.l to prove that 



lim X ~ a, 

T — > + CxD 



where x — (E+,i?, 7i) and a is given by the right-hand sides of (|B.12|I . 
I)B.13|I and ljB.16p . considering the three cases |<7<|,7=| and 
I < 7 ^ 2 simultaneously. 

We begin by defining the subset D in theorem IB. 1 1 bv 

We now verify the hypotheses Hi and H2- Firstly, let w : [tq, +00) ^ D be 
any C^[to,oo) function. Since limT-^+oo M(r) = it follows immediately 
from lprT2|l that 



liin g{w{T),T)= lim M{t) [B^^, RBi,,nB^ 



-0, 

x—w{t) 



showing that Hi is satisfied. Secondly, H2 is satisfied since the variables 
E+, R and H are bounded for all t ^ tq with tq sufficiently large. There- 
fore, since 

lim x{t) = a 

r — *-t-oo 



for all initial conditions a;(To) in D (see (IB.12II . (jB.13|) or (IB.lGp ). theo- 
rem implies that 

lim x{t) = a (B.17) 

r — ^+00 

for all initial conditions x{tq) in D. 

Finally, we need to show that any initial condition a; (tq) — (E+,i?, 7i)|^_^ 
M(to), '0(ro) for the DE ifT^ . subject to > and determines 
an initial condition x{to) in D for the DE H3.11|l . so that ljB.17|) is satisfied. 
Indeed, since limT-^+oo 4' = +00, we can without loss of generality restrict 
the initial condition ■(/'(''o) to be a multiple of tt. This requirement can be 
achieved by simply following the solution determined by the original initial 
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condition until this condition is satisfied. It follows from this condition, in 
conjunction with (|3.8() and the restriction il > applied to (|2.23() . that 

{tl + R' + In') U^^^ ^ (E^ + R' + In') U^^^ < 1, 

so that S(ro) € D. 

Appendix B.3 

We now provide the proof of H3.5|) for the case | < 7 < |, which gives 
the limit of the ratio R/M at late times. In analogy to ()3.8|l . we define a 
variable M by 

M = M (1 + \R'Msm2iP) . (B.18) 
It follows from (|2.21|l that the evolution equation for M is of the form 

M' = -{Q + 2t+ + MBM)M, (B.19) 

where Bjyj is a bounded function for r sufficiently large. By using (|3.9|) we 
obtain _ ^ 

f ) - [3(7 - 1) + M*, M, ^)] (^^) , (B.20) 
where 

h{x, M, ?/') = 2 - 37 + 2Q + 3S+ + MB^ 

and is a bounded function for r sufficiently large. It follows from H3.6|l . 
1)3. 8|1 . H3.10|l and theorcm l3.1l that limT-^+oo h{x,M,ip) — 0. Consequently, 
ljB.20|) implies that 

i? 



^^0(e[3(7-iH^l^), if |^7<1, 
M 

'r 



^ = O (e[3(i-^)+*l-) , if 1<7<|, 



as 



r ^ +c» for any 5 > 0. Therefore, on account of l)B.18(l and 13.8|l . 

R fo, if fs;7<l. 



lim 

T-*+oo M +00, if 1 < 7 < 
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It remains to deduce the limit of R/M as r — > +00 for the case 7=1. 
To proceed we compute the asymptotic form of R and M as t — > +00. 
The calculation parallels that for the non-magnetic Bianchi VIIq models 
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detailed in appendix B of PU" . It follows that any solution of the DE (|2.21|) 
subject to the restrictions (|2.25l) with | < 7 < | satisfies 

7i = CHei/2(37-4)r [i + oie-'n] , 

M=CMel/2(2-37)r [1 + 0(6-"-)] , 

as r — > +00, for some constant & > 0, where Cr, Cn and Cm are positive 
constants which depend on the initial conditions. Therefore, 



as T — > +00 and hence 



lini i^^^^O, if 7=1. 
r— ►+00 M Cm 

Appendix C 

In this appendix we give an expression for the Weyl curvature parameter 
yV in terms of the Hubble-normalized variables E+, R, TC, M and -0. Let 
Eai3 and H^fj be the components of the electric and magnetic parts of the 
Weyl tensor relative to the group invariant frame with Cq = u. It follows 
that Ea_j3 and Hap are diagonal and trace-free and hence they each have 
two independent components. In analogy with we define 

^+ — ^{^22 + £33)1 £- — :r7^{£22 ~ £33), 

1 " 1 (c.i) 

7Y+ = 2(^22 + ^- = 5^(^22-^33), 

where £ap and Ti.af3 are the dimensionless counterparts of Eaf3 and Ha/3, 
defined by 

P Ea(3 ^ Hal3 ,^ „^ 

Oal3 - -Jp^ 'la/3 — 'Jp' ' y^-"^) 

It follows from ((23, ifUTT) and that 

^£l + £l + nl + nl. (C.3) 

Equations (1.101) and (1.102) in for Eap and Hap, in conjunction with 
the frame choice detailed in section 2 of and equations (jC.l|l and (|2.20|) 
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in the present paper lead to 




(C.4) 
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